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The extinction time of an isolated population can be exponentially reduced by a periodic modu- 
lation of its environment. We investigate this effect using, as an example, a stochastic branching- 
annihilation process with a time-dependent branching rate. The population extinction is treated in 
eikonal approximation, where it is described as an instanton trajectory of a proper reaction Hamil- 
tonian. The modulation of the environment perturbs this trajectory and synchronizes it with the 
modulation phase. We calculate the corresponding change in the action along the instanton using 
perturbation techniques supported by numerical calculations. The techniques include a first-order 
theory with respect to the modulation amplitude, a second-order theory in the spirit of the Kapitsa 
pendulum effect, and adiabatic theory valid for low modulation frequencies. 
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I. INTRODUCTION 

Extinction of species after maintaining a long-lived 
quasi-stationary self-regulating population is a striking 
manifestation of large deviations coming from intrinsic 
stochasticity of processes of birth-death type. Not sur- 
prisingly, the extinction phenomenon is at the center 
stage of population biology and epidemiology More 
recently, it has attracted interest in the context of cell 
biochemistry Q. As birth-death processes are usually 
far from equilibrium, they are also of much interest to 
physics Birth-death processes often occur in a 

time- varying environment, and the variations of the en- 
vironment manifest themselves as an explicit time depen- 
dence of the birth and/or death rates. Elucidating dif- 
ferent regimes of exponential reduction of the mean time 
to extinction (MTE) of species due to the time variation 
of the process rates is both important and interesting 
There is a large body of work on the effects of en- 
vironmental noise on the MTE of birth-death systems. 
Early theoretical works on this subject assumed that the 
environmental noise which modulates the rates is delta- 
correlated [1, 0] ■ More recently, effects of finite correla- 
tion time of the noise have also been addressed in nu- 
merical simulations [8[. Not surprisingly, the simulation 
results provide only a partial understanding of the com- 
plex interplay between the nonlinear kinetics and intrin- 
sic stochasticity of the self-regulating population on the 
one hand, and the magnitude and spectral/correlation 
properties of the environmental noise on the other hand. 
To get more insight into this type of problems, we will 
follow the recent work of Escudero and Rodriguez Q and 
consider a much simpler model, where the time variation 
of the environment, as manifested in the rate modula- 
tion, is periodic in time. Although completely ignoring 
the noise aspect, this model enables one to investigate, in 
a controlled and often analytic way, different frequency 
and amplitude regimes of response of a self-regulating 
birth-death system to the rate modulations. Further- 
more, not all environmental variations are noisy: some 
of them (such as daily, monthly, and annual cycles) oc- 



cur in an almost periodic fashion. 

The main objective of this work is to calculate, in dif- 
ferent regimes of parameters, the exponential reduction 
in the MTE of birth-death systems due to a sinusoidal 
rate modulation. These calculations are intimately re- 
lated to finding the "optimal path to extinction" , or in- 
stanton connection in birth-death systems [l(| [H| ■ The 
instanton connection emerges in eikonal approximation 
to the original master equation, and it describes the most 
probable sequence of events that brings the birth-death 
system from its long-lived quasi-stationary state to ex- 
tinction. As a prototypical example we will consider the 
stochastic branching-annihilation process A A 2A and 
2 A 0. This single-species birth-death process can be 
viewed as a simplified version of the well-known Ver- 
hulst logistic model [T3|. Its constant-rate version was 
previously analyzed by a number of methods employ- 
ing the large parameter Q = A//i 3> 1 that corresponds 
to a large average number of individuals in the quasi- 
stationary state at times short compared to the MTE 

\HiMMM- 

The main theoretical tool that we employ in this work 
is an eikonal method which assumes that fl > 1, We 
apply the eikonal method to the exact evolution equa- 
tion for the probability generating function that encodes 
the original continuous-time master equation 0, 0|- In 
the geometrical-optics order of the eikonal approxima- 
tion, the evolution equation for the probability gener- 
ating function reduces to a Hamilton-Jacobi equation, 
with an effective Hamiltonian that explicitly depends on 
time. We explore, analytically and numerically, the phase 
trajectories generated by this Hamiltonian and find the 
perturbed instanton: a perturbed heteroclinic phase tra- 
jectory that corresponds to extinction. Calculating the 
action along the instanton yields an approximation for 
the logarithm of the MTE [l(| EH • Because of the ex- 
plicit time dependence of the eikonal Hamiltonian, an- 
alytical progress is possible only via perturbation ap- 
proaches. Therefore, we use the eikonal method in con- 
junction with two different perturbation techniques that 
employ additional small parameters. The first of them 
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is linear theory (LT): a theory based on linearization of 
Hamilton equations with respect to the modulation am- 
plitude. For birth-death processes, the LT has been re- 
cently used by Escudero and Rodriguez A closely 
related theory was earlier developed by Dykman et at 
(l6l in the context of large fluctuations and escape 
in continuous (Langevin-type) stochastic systems driven 
by a time-periodic signal. Dykman et al. showed that in 
LT the modulation signal removes the degeneracy of the 
unperturbed instanton trajectories with respect to time 
shift. This leads to the selection of the optimal instanton 
which is synchronized with the phase of the modulation 
signal. In general, the action along the optimal instanton 
is smaller than that the action along the unperturbed in- 
stanton, leading to an exponential reduction of the escape 
time (Mini- 

By employing a theorem due to Melnikov QUEi, Es- 
cudero and Rodriguez [9( proved the existence of a per- 
turbed instanton in the branching-annihilation process 
with weakly modulated reaction rates. We extend the 
LT by calculating analytically the action along the in- 
stanton, which yields the logarithm of MTE. We also 
show numerically that the perturbed instanton persists 
for any reasonable modulation amplitudes, and for all 
modulation frequencies. Furthermore, we show that the 
LT gives satisfactory results only when the modulation 
frequency u> is smaller than or comparable with the relax- 
ation rate of the system A. When lu 3> A, the first-order 
correction to action, as predicted by the LT, turns out to 
be exponentially small with respect to the rescaled fre- 
quency lu/X. It is the second-order correction in the mod- 
ulation amplitude that becomes dominant in this regime. 
We calculate this correction by employing (a Hamiltonian 
extension of) the Kapitsa method, see e.g. [2fj| . 

The opposite, low-frequency regime is both the sim- 
plest and most important as, for the same modulation 
amplitude, the exponential reduction of the MTE turns 
out to be the largest here. This regime can be efficiently 
dealt with in adiabatic approximation, for any reason- 
able modulation amplitude. Here one assumes that the 
extinction rate corresponding to the unperturbed prob- 
lem is known, obtains the instantaneous extinction rate 
and averages it over the modulation period. This simple 
procedure yields the average extinction rate, and hence 
the MTE, of the perturbed problem. If the knowledge 
of the extinction rate of the unperturbed problem in- 
cudes a pre-exponent, the adiabatic approximation yields 
a (modified) pre-exponent of the MTE of the perturbed 
problem: a significant improvement over the geometrical- 
optics order of the eikonal method. 

The remainder of the paper is organized as fol- 
lows. We begin Section |TT] with a brief introduction 
to the geometrical-optics theory of extinction in birth- 
death processes with time-independent rates using the 
branching-annihilation model as an example. In Sec- 
tion IIIII we present the eikonal theory of population ex- 
tinction for a time-periodic rate modulation. We begin 
Section IIIII by developing the LT, and continue by ad- 



dressing the high-frequency limit lu 3> A and employing 
the Kapitsa method. We conclude Section [TIT] by report- 
ing numerical solutions of the Hamilton equations, in or- 
der to verify our theoretical results and extend them be- 
yond the validity domain of the perturbation techniques. 
Section lTVl deals with the low-frequency limit by means of 
adiabatic approximation. Theoretical results obtained in 
this limit are compared with numerical solutions of the 
master equation. Section [V] presents a brief summary 
and discussion of our results. 



II. MASTER EQUATION, PROBABILITY 

GENERATING FUNCTION, AND 
GEOMETRICAL OPTICS OF EXTINCTION 

Here we present a brief introduction into an eikonal 
theory of population extinction in a time-independent 
environment [lOL fllj , using the prototypical example of 
the branching-annihilation process A A-s- 2 A and 2 A ^ 0, 
where A, (1 > are the reaction rate constants [ill EE 
Hil . 0j|. We assume in this section that A and \x are 
independent of time. 

At the level of deterministic modeling, the dynamics 
of the average number of individuals n(t) is described by 
the (mean-field) rate equation 

_ = An- M n. (1) 

This equation predicts an attracting fixed point n = 
X/fi = fl 1. The rate equation ignores the intrin- 
sic noise that comes from the discreteness of individuals 
and stochastic character of the reactions. This noise de- 
termines the probability distribution of the actual values 
of n and, in particular, a quasi-stationary distribution 
around the average value f2 3> 1 that sets in on a relax- 
ation time scale A -1 , see below. On this time scale the 
mean- field picture correctly describes the average number 
of individuals in the stochastic process. At sufficiently 
long times, however, the noise invalidates the mean-field 
picture completely due to the existence, in the stochastic 
process, of the absorbing state n = 0. In other words, 
the stochastic process A — > 2A and 2A — > eventually 
suffers a rare sequence of events that brings the system 
into the empty state. 

The intrinsic noise is accounted for quantitatively by 
the master equation which describes the time evolution of 
V n {t): the probability to have n individuals at time t. In 
the branching-annihilation example, the continuous-time 
master equation, for n > 1, reads 

= !t{(n + 2)(n + l)V n+2 (t)-n(n-l)V n (t)} 
+ A [(n - l)V n -i(t) - nV n (t)] . (2) 

The master equation can be conveniently recast with the 
help of the probability generating function 0, 0] 

00 

G(p,t)=Y / P n 'Pn(t), (3) 
n=0 
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where p is an auxiliary variable. G(p, t) encodes all the 
probabilities V n (t), as the latter ones are given by the 
coefficients of the Taylor expansion of G(p, t) around 
p = 0. As the probabilities are normalizable to 1, 
Y^o''Pn(t) = 1, the generating function satisfies the con- 
dition 



G(l,t) = 1. 



(4) 



The time-dependent moments of the distribution V n (t) 
can be expressed through the derivatives of the gener- 
ating function at p = 1, e.g. (n)(t) = *^2 n nV n (t) = 

After a simple algebra, Eq. ^ can be transformed into 
an exact evolution equation for G{p,t): 

dG n 2 .8 2 G . . dG 

The absorbing state at n — corresponds to the station- 
ary solution G s (p) = 1. The presence of the absorbing 
state is manifested by the absence in Eq. ([5]) of a term 
proportional to G(p, t). 

An initial value problem for Eq. §5§ can be solved by 
expanding G(p, t) in the eigenmodes of a Sturm-Liouville 
problem related to the non-Hermitian operator 



d 



dp 



dp 



see Ref. [15[ for detail. As in many other problems of 
extinction of species, there is a zero eigenvalue Eq = 
which describes the absorbing state, and an infinite dis- 
crete set of eigenvalues {E n }%L 1 which describe an expo- 
nential decay with time of the rest of eigenmodes (l5l.[2l|. 
For fl ^> 1 the eigenvalue Ei is exponentially small in 
O, and much less than the eigenvalues E 2 , E%, . . . , which 
are of the order of the relaxation rate of the system A. 
Therefore, there are two widely different time scales in 
the problem. The decay, on the time scale A -1 , of the 
higher eigenmodes corresponds to a rapid relaxation of 
the probability distribution V n (t) to the quasi-stationary 
distribution, peaked at n ~ f2. The much slower decay, 
on the time scale E^ , of the lowest excited eigenmode 
[accompanied by a slow growth of the extinction proba- 
bility Vo(t)} corresponds to an exponentially slow decay 
of the quasi-stationary distribution and to extinction of 
the species. A simple approximation to the eigenvalue 
Ei is provided by geometrical optics, by which we mean 
the leading order of the eikonal approximation [ill ] . Al- 
though it gives only exponential accuracy, the geometric 
optics can be readily applied to systems of two species 
[22l [23l ] and to time-dependent problems [H, [H, [2J] : the 
focus of the present work. 

Employing the eikonal ansatz G(p,t) = exp[— S(p, t)} 
in Eq. (O and neglecting d 2 S/dp 2 , we arrive at 
a Hamilton-Jacobi 
representation: 



equation for S{p,t) in the p- 



dS ii, 2 , (dSy x , ,dS „ 

m + l il -^{d-p) -ap(p-D^ = o. (6) 




FIG. 1: (Color online.) The phase plane of the unperturbed 
Hamiltonian ((7}. The solid lines denote zero-energy trajecto- 
ries. 



Introducing a canonically conjugate coordinate q — 
—OS/ dp and shifting the momentum p — p— 1, we arrive 
at a one-dimensional Hamiltonian flow, where p plays the 
role of the momentum [ill: 



H(q,p)= X(l+p)-^(2+p)q 



qp. 



The Hamilton equations are 
OH 

Q = -g- = q[X(l + 2p)-fi(l+p)q], 
P = —5- =p[»(2+p)q- A(l+p)] 



(7) 



(8) 
(9) 



As the geometrical-optics Hamiltonian ([7]) does not de- 
pend explicitly on time, it is conserved: H{q,p) = E = 
const, and the problem is integrable. The phase plane 
(q,p), defined by this Hamiltonian flow, provides a useful 
visualization of the dynamics, see Fig. [TJ The attracting 
fixed point n — Q and the repelling fixed point h = 
of the one-dimensional rate equation ((TJ) become hyper- 
bolic points A(q = n, p = 0) and 0(q = 0,p = 0) of the 
phase plane (q,p). They belong to the mean- field zero- 
energy line p = 0. [The equality H(q 7 0) = reflects the 
probability conservation.] Another zero-energy line is the 
extinction line q = 0. Of special importance, however, is 
the non-trivial zero-energy line 



q = qo(p) = 



(2+P) 



(10) 



This line, which includes still another hyperbolic fixed 
point B(q = 0, p = —1) (we call it the fluctua- 
tional point), gives a geometrical-optics description of the 
quasi-stationary probability distribution. The — 1 < p < 
segment of this line is a heteroclinic trajectory. It exits, 
at time t = — 00, the hyperbolic fixed point A along its 
unstable manifold and enters, at time t = 00, the hyper- 
bolic fixed point B(q = 0, p = —1) along its stable mani- 
fold. The heteroclinic trajectory (often referred to as the 
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optimal path to extinction, or the instanton [25l. l2a. |27| ) 
describes the most probable sequence of events bringing 
the system from the quasi-stationary state to extinction. 
Up to a pre-exponent, that is missed by the geometrical- 
optics order of the eikonal theory, the MTE can be ap- 
proximated as r ex ~ exp(£>o), where 



S 



q (p) dp = 20(1 - In 2) 



The time-dependent instanton coordinate qo(t 
momentum po(t — to) are the following: 



(11) 



to) and 



<?o (t - t ) = 



2fi 



2 _|_ e X(t-t ) 



Po(t - t ) = -- 



1 + e- A (*-*o) ' 
(12) 

where to is an arbitrary time shift. 

The eikonal approximation is valid when Sq > 1, 
that is f2 3> 1. Equation (fITj) . obtained by Elgart and 
Kamenev [ll|, coincides with the leading order term of 
the logarithm of the MTE found by other methods (also 
employing the strong inequality 1! > 1), by Turner and 
Malek-Manso ur J13| . Kessler and Shnerb [14]], and Assaf 
and Meerson [la ]. 



III. EXTINCTION IN A TIME- VARYING 
ENVIRONMENT: PERTURBATION 
TECHNIQUES FOR THE INSTANTON 

Now let the branching rate A in the master equation 
([2]), and in the evolution equation ([5]) for G(p,t), vary 
sinusoidally in time: 



A(i) = A (1 + ecosujt) , 



(13) 



where \e\ < 1, and e and u> are the modulation amplitude 
and frequency, respectively. The eikonal approximation 
now yields the following time- dependent Hamiltonian: 



where 



and 



H(q,p,t)=H Q (q,p)+eH 1 (q,p,t), (14) 



H (q,p) = A (l +p)pq - ^-(2+p)pq 2 , (15) 



Hi (q, p, t) = A (l + p)pq cos ut . 



(16) 



For the unperturbed Hamiltonian Ho(q,p), the instanton 
(fTT)|) and the corresponding action So along it, Eq. (fTTjl . 
are known. Our strategy will be to calculate the action 
along the perturbed instanton: the instanton of the per- 
turbed, time-dependent Hamiltonian (Ti"4"]) . 

Following Escudero and Rodriguez [9(, let us consider 
the Poincare map P*° = {q,p,t\t = t S [0,T]}: the 
projection of three-dimensional trajectories of the non- 
autonomous system (|14p in the (q,p, t) space on the (q,p) 



plane at section time to € [0, 2ir/tu} . In the perturbed sys- 
tem the hyperbolic points A and B are also perturbed. 
We will denote the perturbed fixed points by A^ and 
B*° , respectively. For generic Hamiltonians the existence 
of the perturbed hyperbolic points is guar anteed, by the 
Poincare-Birkhoff fixed point theorem , only for suffi- 
ciently small e. For the birth-death Hamiltonian we are 
dealing with here the perturbed fixed points exist for any 
\e\ < 1. The point -B' = (q = Q,p = —1) remains un- 
changed [28[, whereas the point A l ° can be easily found 
from the time-dependent rate equation ([5]) with p = 
and A from Eq. (Ill: 



(17) 



q = q [Aq ( 1 + £ cos ut) — /iq] 



The solution of this equation, for an arbitrary initial con- 
dition q(0), is the following: 



?(*) 



q(0) exp (\ t 



eXn sin ujt \ 



^(0)/ o Aot exp 



s + 



eXo sin(ujs/Ao) 



ds 



, (18) 



see also Ref. @. At long times this solution becomes 
periodic in time, and can be written as 



lit) 



exp ( 



eXp sin ut 



/-^O I n (^ a ) COs(nLUt — <j> n — U7T /2) 



, (19) 



n— — oo 



where /„(...) is the modified Bessel function, and <p n = 
arctan(nw/Ao). Putting in Eq. (| 1 9)) t = to we obtain q 
of the perturbed fixed point A' , whereas its p remains 
zero. 

Existence of the perturbed fixed points is a necessary 
but, in general, insufficient condition for the existence of 
a hetcroclinic trajectory connecting the unstable mani- 
fold of A* and the stable manifold of £?*° . In general, one 
should also establish the existence, in the perturbed prob- 
lem, of the unstable and stable manifolds themselves, 
and of their intersection 0, . We verified analytically 
and numerically, see below, that the unstable manifold 
of A l ° and the stable manifold of £?*° do intersect. Let 
us denote the perturbed instanton connection by the pair 
q(t, to),p(t, to). As the energy is not conserved, the action 
along the perturbed instanton now includes an integral 
of H over time: 



S = / {p(t,t )q(t,to)-H [q(t,to),p{t,t )} 
- eJri[g(t,to),p(*,to),t]}dt, (20) 

where q{t, to) = dq/dt, and the unperturbed Hamiltonian 
Ho{q,p) is invariant to the specific choice of to- Escudero 
and Rodriguez 0] applied the Melnikov theorem [3, [l9[ 
and proved that a perturbed instanton exists in this prob- 
lem for sufficiently small e. To this end they calculated 
(an approximation for) the distance between the unsta- 
ble and stable manifolds, which is given by the Melnikov 
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function, see below, and showed that it vanishes for a spe- 
cific choice of to . We briefly reproduce their derivation in 
subsection IIII Al along with new results: a linear-theory 
calculation of the action (j2"0|) and the corresponding re- 
duction in the logarithm of the MTE. We also show, in 
subsection IIII Bl that the perturbed instanton exists for 
high-frequency perturbations, w > Ao at any |e| < 1. 
Furthermore, we report, in subsection IIII Cl strong nu- 
merical evidence that the perturbed instanton exists for 
any |e| < 1 and any lu. 



A. Linear correction to action 

Here we consider the linear theory (LT) which assumes 
that the term eH\ (q, p, t) in Eq. (fTij) can be treated per- 
turbatively. For this assumption to hold it is sufficient to 
demand the strong inequality \e\ <C 1. As we will see in 
subsection IIII B | at high frequencies, ui ^> Xq, the strong 
inequality \e\ -C 1 becomes unnecessary. For complete- 
ness, we will derive Eqs. (l24l) and (f25|) in a general form, 
before dealing with the particular example of branching 
and annihilation. In the first order in e the perturbed 
instanton of H(q,p,t) is described by the equations 

q(t,t ) = £fo(t-to) +eqi(t,t ), 

p(t,t ) = po(t-<o) + epi(t,*o), (21) 

where qo(t — to) and po(t — to) stand for the (known) 
instanton solution of the unperturbed equations 



<?o = 
Po = 



dHo(q Q ,p ) 

dpo 
dHo(qo,po) 
dqo 



(22) 



that is Eqs. © and © with A = A . 

To calculate the action (f2T)|) we expand the integrand 
in e and obtain, in the first order, 



(Po +epi)(q + eqi) - H Q (q ,po) - eq x 



dH (q ,Po) 



dH (q ,Po) rj , ,s ■ tt i \ 

epi ~ eHi{qo,p ,t) ~ p qo - H (q ,Po) 

opo 

+ spiQo + £Poqi + £QiPo ~ spiQo ~ e^i(9o,Po, t) . (23) 

After the integration the first two terms yield the un- 
perturbed action So- The third term cancels out with 
the sixth term, while the fourth term cancels out, after 
integration by parts, with the fifth term. The result is 
S(to) = So + AS (to), where 



AS (t Q 



#i[?o(i-io),Po(i-*o),*] dt. (24) 



Recall that the integration is performed along the un- 
perturbed instanton. To find the optimal correction to 



action we have to minimize 5 (to) with respect to to (com- 
pare to Refs. [H,[l3])- This yields the following equation 
for to: 



dS(t ) 
alto 



= e 



— oo 
oo 



~5 — °a 



po dt 



—e 



dpo 

/OO 
{H ,Hi} dt = 0, 
-OO 



(25) 



where {Hq,Hi}o is the Poisson bracket evaluated on 
the unperturbed instanton. The quantity M(to) = 
{Ho,Hi} Q dt is the Melnikov function of the per- 
turbed problem [3 . It is proportional to the dis- 
tance between the unstable and stable manifolds of Al° 
and -B* , respectively [lj|. That M(to) has simple zeros 
yields a sufficient condition for the existence of the per- 
turbed instanton. These zeros are the critical points of 
the function S (to)- By finding the critical value of to for 
which S (to) has its minimum, we obtain the minimal ac- 
tion along the instanton of the perturbed problem (|24|) . 
It is the min AS (to) that yields the reduction of the log- 
arithm of the MTE, caused by the time-dependence of 
the rate. Note that when Ho commutes with Hi, M(to) 
vanishes identically at any to , but this degenerate case is 
hardly of interest. 

The validity of the linear eikonal theory is determined 
by the double strong inequality 



1 <C AS -C So 



(26) 



The left inequality is required for the eikonal approxima- 
tion to be valid, while the right inequality is needed for 
the linearization to be valid. Note that the rate variation 
does not have to be time-periodic for the LT to hold: it 
suffices to demand that the rate variation be bounded by 
a sufficiently small value. 

Now we return to the branching-annihilation example. 
In the first order in e, the perturbed fixed point Al° is 
determined by the equation 



= 1 



e\o (Xo cos tot + uj sin ujt) 

^Tx* 



(27) 



evaluated at t = to- Substituting Eqs. (fT2|) and (fl6|) in 
Eq. pgj). we obtain 

Calculating the integral we arrive at 



S(t ) 
So 



x < — COSLUto 

Xo 



(1 - In2)sinh(7rw/A ) 
'In 2 



sin 



A 



■+*t 



sin uit 



(29) 



where So is given by Eq. The minimal action is 

S(t* Q ), where tp is the solution of the trigonometric equa- 
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FIG. 2: (Color online.) The difference uta — n, see Eq. pi) , 
as a function of a = lo/Xq. 



tion for to'. 

^ - -eM(t ) = 

dtn 



x < — sin iota + cos 
Xq 



2ir eQlu 

sinh(7njj/Ao) 

In 2 
w I — hto 



— cos ojto } = 0, 
(30) 



subject to the condition d 2 S(to) / 'dt^ > 0. The solution 
is 



cjf = 7r + arctan 



1 — cos(aln2) 



a — sin(a In 2) 



(31) 



where a — lu/ Xq, and we assume here and in the remain- 
der of this subsection that e > [3(j. Equation ([3"T]) 
shows that the optimal instanton becomes synchronized 
with the rate modulation, similarly to what happens in 
the context of escape in Langevin-type stochastic systems 
driven by a time-periodic signal [Ig,[l7|]- Figure[5]depicts 
the difference tuto — 7r as a function of a. This difference 
is small at a < 1, u>t - it ~ (ln2) 2 /[2(l - ln2)]a, and 
at a > 1: ujt — ir ~ a _1 [l — cos(aln2)]. ft vanishes 
for a = 27rn/ln2, where n = 1,2,.... For a ~ 1 the 
difference w^o — T is of the order of one. 

Now we use Eqs. (|29|) and (|3Tj) to find the perturbed 
action: 



5 



1 



AS 



1 



7T£ 



(1 - In 2) sinh(Tra) 
|[sin(aln2) - a} 2 + [cos(aln2) - 1] 2 | 



1/2 



(32) 



Figure [3] depicts AS/ (efi) as a function of a as predicted 
by Eq. (|3"2"]) . The maximum effect of the rate modulation 
is obtained at a — > 0. Recall that the perturbed action 
yields (in the leading order in f2) the natural logarithm of 
the MTE of the perturbed branching-annihilation prob- 
lem. 
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FIG. 3: (Color online.) The correction to action A<5>, rescaled 
by eSo, versus the rescaled modulation frequency a = lo/Xq. 
Dashed and dash-dotted lines: predictions of the linear the- 
ory, ASlt, see Eq. (|32[l . for e = 0.04 and e = 0.1, respectively. 
Solid line: AS num found by solving the Hamilton equations 
numerically, sec subscction lHI CI On this scale the three lines 
are indistinguishable. The inset shows the ratios of ASlt and 
ASnum for e = 0.04 (the dashed line), and e = 0.1 (the dash- 
dotted line). The disagreement at high frequencies demands 
an account of what we call the Kapitsa correction, calculated 
in Section IXTTBl 



The perturbative eikonal result is valid as long as the 
double inequality (|2l)|) is obeyed. The left inequality 
breaks down at very small e, whereas the right inequal- 
ity breaks down at not small e where linearization be- 
comes invalid. One can define, in analogy with the prob- 
lem of activated escape in Langevin-type stochastic sys- 
tems the logarithmic susceptibility x = d(AS)/d\e\. 
Within the framework of the LT [that is, for intermedi- 
ate values of the modulation amplitude that satisfy the 
double inequality (|26p]. x is independent of e: 



X = 



sinh(7ra) 

|[sin(aln2)-a] 2 + [cos(aln2)-l] 2 . (33) 



At small and large modulation amplitudes x becomes 
amplitude-dependent. 

In the low- frequency limit, a <C 1, one obtains from 
Eq. ~ 



AS a 2 
eSo 6 



7T 2 - 



(ln2) 



(ln2) 4 



(l-ln2) 2 4(1 -In 2) 



, (34) 



where we have kept the leading and subleading terms. 
The leading term coincides with the prediction ([55]) from 
the adiabatic approximation, see Sect ion HVl In this limit 
the LT is valid for il^ 1 <C e <C 1. The same criterion 
applies for a ~ 1, as AS/(eVt) in Eq. (|3"2"|) is of order 1 
there. In the high-frequency limit, a»l, the correction 
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AS is exponentially small: 



order in £ and 77: 



AS 

eS 



2na e 



l-ln2 



(35) 



Here e need not be small compared to 1 to satisfy the 
right inequality in Eq. {26]). Note that Eqs. p2 ]) - ([35]) 
with e replaced by |e| are valid also for e < 0. 

In Fig. [3] Eq. (f3"2"|) is compared with the correction to 
action found by solving the Hamilton equations numer- 
ically (see subsection IIII CI for details of the numerics). 
One can see a good agreement at small and intermedi- 
ate frequencies, and a disagreement at high frequencies. 
What is the reason for the disagreement? That the cor- 
rection to action AS, as predicted by the LT, becomes 
exponentially small at large frequencies indicates that the 
e 2 -correction, neglected by the LT, becomes dominant 
there. This conjecture is validated in the next subsection, 
where we consider the high-frequency case and calculate 
the £ 2 -correction using (a Hamiltonian extension of) the 
Kapitsa method. 



H(t , P , t ^H(x,Y, )+i a Jim + v 9H(x - Y) 



dX 



dY 



£ 2 d 2 H(X,Y) t d 2 H(X,Y) + ^d 2 H{X,Y) 



2 dX 2 
ee H(X,Y,t). 



8Y 2 



dXdY 



(37) 



The Hamilton equations become 
dH(X,Y,t) 



q = X + £ 



dY 



p = Y + 77 ~ — 



dH(X,Y,t) 
dX 



(38) 

Now we demand that the rapidly oscillating terms in 
Eqs. (|38|) balance each other. This yields 



£ ~ e\ X(2Y + 1) cos wt , 7) ~ -eX Y(Y + 1) cosut , 

where terms of the order of £ and 77 have been neglected, 
but their time derivatives (which are proportional to ui 
and therefore large) have been kept. Treating X and Y 
as constants during the period of rapid oscillations 2ir /to, 
we can easily solve Eqs. (f3U)) and obtain 



B. Quadratic correction to action: the Kapitsa 
correction 



Here we consider the high-frequency limit, to Xq 
and calculate the correction to action which is of the 
second order in e. To this end we use (a Hamiltonian 
extension of) the method that was developed long ago 
in the context of the "Kapitsa pendulum" : a rigid pen- 
dulum with a rapidly vibrating pivot, see e.g. [20]. The 
Kapitsa method involves, as a first step, calculation of 
small high-frequency corrections to the unperturbed co- 
ordinate and momentum of the system. Because of the 
high frequency of the perturbation these corrections are 
small even if e is of order 1. Using the high-frequency cor- 
rections, we construct a canonical transformation which, 
by means of time averaging (that is, by rectifying the 
high-frequency component of the motion), transforms 
the original time-dependent Hamiltonian into an effec- 
tive time-independent one. The effective Hamiltonian 
includes a correction coming from the rectified high- 
frequency perturbation. Finally, we find the perturbed 
instanton, emerging from this effective time-independent 
Hamiltonian, and the action along the instanton. 

The starting point of the derivation is the same Hamil- 
tonian ([21) with rates given by Eq. (|13p . We represent q 
and p as follows: 



q(t)=X(t)+t(t), p(t) = Y(t) + r,(t), 



(36) 



where X and Y are slow variables, while £ and r\ are small 
and rapidly oscillating. Now we expand the Hamiltonian 
H(q,p,t) around q = X and p = Y up to the second 



eX 



X(2Y+l)smut , rj(t) 



eX 



r(y+l)sinwt. 

LU LU 

( 40 ) 

Now it is clear that this perturbation scheme demands 
|e|Ao/w <C 1. As we have assumed u> 3> Ao, £ need not 
be small. 

Using Eqs. (|36| and (|40|) . we perform an almost canon- 
ical transformation from the old variables q and p to the 
new variables X and Y: 

q(X,Y,t) 



1 _ £^o(2F + I) smut 



fX 

1 + — {2Y+l)smujt- 

LU 



~ X 
p(X,Y,t)=Y 



^°(2y + l) 2 sin 2 ^ 



1 _ £^(i + y) sincji 

LU 



(41) 



This transformation is canonical up to third order of 
0[(Xq/lu) 3 ] <C 1, as the Poisson brackets {q,p}(x,Y) 
1 + C[(Ao /w) 3 ] ■ The generating function of this transfor- 
mation is 

F 2 (q, Y, t) =qY[l- {eX /uj){Y + 1) sin cut] , 

see e.g. Ref. [20| . Now we transform to the new variables 
X and Y, H' = H + dF 2 /dt, and average the new Hamil- 
tonian H 1 over the period of rapid oscillations 2tt/uj. This 
yields an effective Hamiltonian 

H(X, Y) = X n XY{l + Y) - — (Y + 2)YX 2 



F 2 \ 2 



X Y(l + Y) 2 



/'o 



(2 + 12Y + 18Y 2 + 5Y 3 ) 



(42) 
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which is time-independent. The first two terms come 
from the unperturbed Hamiltonian, Eq. (|15p . the other 
two terms describe a rectified e 2 -correction coming from 
the high-frequency perturbation. The mean-field fixed 
point and fluctuational fixed point of the effective Hamil- 
tonian are 



A 



Q 1 



2u 2 



p = 



and B(q = 0, p = —1) , 



respectively. As the effective Hamiltonian is time- 
independent, there is no need in the Poincare section in 
this limit. Two "trivial" zero-energy lines of the effective 
Hamiltonian are X = and Y = 0. The non-trivial zero- 
energy line yields the effective instanton Xq(Y). In the 
leading order of a = w/Xq ^> 1 we obtain 



Xo(Y) 



2fi(i + y) 

2 + Y 



1 - 



e 2 (2 + ior + i5y 2 + 4r 3 ) 



2a 2 (2 + Y) 



The action along the effective instanton is 



S = - Xo(Y)dY = S (l- 



where 



R ^ 6 In 2 -49/12 
1 - In 2 



(43) 



(44) 



(45) 



We will call the resulting correction to the unperturbed 
action. 



AS K = -KS (e/a) 2 , 



(46) 



the Kapitsa correction. The eikonal approximation, that 
has led to Eq. (|4*6")) . demands ASk S> 1. Figure H] shows 
the ratio of ASk and the correction to action, found 
by solving the Hamilton equations numerically (see sub- 
section UTTC] for details of the numerics). As expected, 
good agreement is observed at high frequencies, where 
the Kapitsa method is applicable. 

Now we can put the results obtained with the LT and 
with the Kapitsa method into a broader context. The 
correction to the unperturbed action AS = S — So < 
must have the following general form: 



AS 



-h{a)\e\-h{a)e 2 + 



(47) 



(Clearly, AS must be analytical in |e|, rather than in e, 
as changing the sign of e only brings about a phase shift 
7r which changes the minimum point to but leaves the 
minimum action unchanged.) The function f\(a) > is 
given by the LT, see Eq. 



A (a) 



n | [sin(a In 2) - a] 2 + [cos(aln2) - 1] 2 | 



(1 - In 2) sinh(Tra) 



1/2 



(48) 



0.75 



CO 

< 



CO 

< 




0.25 



FIG. 4: (Color online.) Fhe ratio of the Kapitsa correction 
AS K , Eq. ((ill, and AS num = 5„ um -2n(l-ln2) for e = 0.1 
(dashed line) and e = 0.2 (dash-dotted line), versus a — lo/Xq. 
The agreement improves with an increase of a and e. 



The high-frequency asymptote of /2(a) is given by 
Eq. 105): 



/2(a) 



K 



a > 1 



(49) 



To calculate the function /2(a) analytically for all a 
would be quite cumbersome, so we found it numerically 
by solving the Hamilton equations for three different val- 
ues of e and for different frequencies (see subsection IIII CI 
for details of the numerics). For each set of e and a 
we computed AS and used Eq. (|47|) with f\{a) from 
Eq. |48|) to extract /2(a)- As expected, the resulting 
plots of /2(a) for different e collapse into a single curve, 
see Fig. [SJ Noticeable is a non-monotonic, alternating- 
sign a-dependence. The high-frequency asymptote of 
/2(a) is in excellent agreement with the Kapitsa correc- 
tion (|49l) . Note that /2(a) vanishes as a — > 0. The reason 
is that, as lu — > 0, the dependence of the eikonal part of 
ASq on e becomes linear, see Section HVl This is a non- 
generic property of the case when the branching rate A 
is modulated, while the annihilation rate /1 is not. 

When does the Kapitsa correction to action dominate 
over the exponentially small first-order correction at high 
frequencies? The corresponding condition e 2 f2(a) S> 
|e|/i(a) can be rewritten as |e| 3> a 3 e _7rQ ; it is easily 
satisfied for large a. In this regime the logarithmic sus- 
ceptibility x is proportional to |e|. 



Numerical calculations of the perturbed 
instanton and the action 



To verify our analytical results, and to explore the pa- 
rameter regions beyond the validity of the perturbation 
methods, we computed the action numerically by solv- 
ing the Hamilton equations. Rescaling time r = Xgt and 
the coordinate Q = /loq/Xo — q/fl in Eqs. (fT^|l - p^|) . one 
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a 



FIG. 5: (Color online.) Shown is the function /2(a) deter- 
mined numerically, see text. The curves for e = 0.04 (solid 
line), e = 0.1 (dashed line), and e = 0.2 (dash-dotted line) 
collapse into a single curve as expected from Eq. (I47p . The 
a 3> 1 asymptote of /2(a), given by Eq. (|49[1 (thick solid 
line), agrees with the numerical result at high frequencies. 
The inset shows a blowup of the high-frequency region. 

obtains rescaled Hamilton equations 
go 

= (l + 2p)Q- (1+ p)Q 2 + e{l + 2 P )Q cos or, 

OT 

dp 

tt- = -(l+p)p + (2 + p)pQ - e(l+ p)p cos ar ,(50) 
ar 

which follow from the rescaled Hamiltonian 

h{Q,p,r) = {p + l)pQ - ^(2 +p)pQ 2 + e(p+l)pQ cos ar. 

The action along the perturbed instanton [Q*{t),p*(t)] is 
given by 

{-Q*(r)p*(r) - h[Q*(T),p*{r),T]}dT. 

-OO 

where the integrand is fully determined by two dimen- 
sionless parameters: the rescaled modulation amplitude 
e and rescaled frequency a — lu/X. The numerical solu- 
tion is obtained by shooting. We start at r = at a point 
[Q(0),p(0)] lying on a small circle of radius S -C 1 cen- 
tered in the unperturbed mean- field point A(Q — l,p = 
0). The polar angle 9 of the point [Q(0),p(0)] on the 
circle serves as the only shooting parameter: it is var- 
ied until the numerical trajectory reaches, at some time 
t = Tf, (a close proximity of) the fluctuational point 
B(Q = — l,p = 0). For the numerically found instan- 
ton [Q*(t),p*(t)] to be a good approximation to the true 
perturbed instanton in its entirety, S must be sufficiently 
small. On the other hand, too a small 5 causes a long r/ 
(and, therefore, a long computation time and possible ac- 
cumulation of numerical errors) because of the intrinsic 
logarithmic slowdown near the fixed point. By repeat- 
ing the computations for a smaller circle, 5\ < S, we 
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FIG. 6: (Color online.) A non-perturbative instanton trajec- 
tory obtained by solving numerically the rescaled Hamilton 
equations (|50p for e = 0.6 and a = 1. The trajectory first 
performs large-amplitude oscillations around the mean-field 
fixed point A of the unperturbed system [or, in other words, 
stays close to the perturbed fixed point of the Poincare map 
A s (to), then leaves the vicinity of the fixed point and ulti- 
mately reaches (a close vicinity of) the fluctuational point 
B(Q = 0, p = —1) which, for the branching-annihilation re- 
action, remains unperturbed. 

checked that the same numerical instanton trajectory is 
reconstructed with high accuracy, apart from additional 
oscillations that appear in the vicinity of the mean-field 
point of the unperturbed system (see Figs. [5] and [7]). As 
long as 5 is chosen to be sufficiently small, the contri- 
bution of these oscillations to the net action along the 
instanton is small. 

We used this algorithm to verify our perturbative the- 
oretical results, see Figs. [3J5] We also computed the 
perturbed instanton numerically in the parameter region 
beyond the validity of the perturbation techniques: e < 1 
and a < 1. These computations strongly suggest that 
the instanton connection persists for any s < 1 and any 
a. Figures [S] and [7] show examples of numerically found 
instantons, and action along them, for non-perturbative 
values of parameters. 

IV. ADIABATIC APPROXIMATION 

Adiabatic approximation is non-perturbative in the 
modulation amplitude e, and demands that the modu- 
lation frequency be much smaller than the characteristic 
relaxation rate of the system. Consider first a general 
single-species birth-death process with time-independent 
rates which possesses a non-trivial quasi-stationary state 
and a single absorbing state at zero. At times much 
longer than the time of relaxation to the quasi-stationary 
state, the extinction probability at time t can be repre- 
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FIG. 7: (Color online.) Projections on the (Q,p)-plane of 
instantons found numerically for four sets of parameters (the 
dashed lines). The parameters and the corresponding rescaled 
actions [see Eq. (|51fl ] S jVL are: (a) e = 0.6, a = 1, S/Q, = 
0.467 (b) e = 0.9, a = 1, <S/fi = 0.398 (c) e = 0.6, a = 
2, S/fi = 0.576, and (d) e = 0.9, a = 2, S/fi = 0.545. 
The solid line denotes the unperturbed instanton (110[) . For 
comparison, the unperturbed rescaled action is <5>o/f2 = 2(1 — 
In 2) = 0.6137 ... . The blowups in the vicinity of p = (the 
insets) show that Q(p) oscillates before departing toward the 
fluctuational fixed point. In these examples S — 10 -6 . 



sented as 



V (t) = \-Ce 



where r® x = r^(X%, A2, ■ ■ • ) is the extinction rate: the 
lowest excited eigenvalue E\, see e.g. Refs. [HI, [2l| . 
The extinction rate depends on the reaction rates of the 
problem Ai, A2, . . . . The constant C is determined by the 
initial conditions; for a sufficientl y la rge initial number of 
particles, C is close to unity (l5ll2l|. 

Now we introduce an adiabatically slow variation into 
the reaction rates so that the characteristic time of the 
variation is much longer than the relaxation time of the 
system, but much shorter than the MTE. The extinction 
probability can now be written as 



P (t) = 1-Ce-ti r -(*')dt' 



(52) 



where r ex (t) = r ex [\i(t), A2(t), • • •] is the instantaneous 
value of the slowly time-dependent extinction rate. The 
average extinction rate f ex can be defined via the relation 



P (t) = 1 - Ce- f ^ . 
Comparing Eqs. (|52")l and (|55|) . we obtain 



= ^J r ex (t')dt', 



(53) 



(54) 



where T is much longer than the relaxation time but 
much shorter than the MTE. For slow periodic rate mod- 
ulations with frequency u Eq. (|54p can be rewritten as 



(55) 



2tt r^/w 
r ex = — / r ex (t') dt' . 
u Jo 



The MTE is equal, in the adiabatic approximation, to 

l/fex- 

Now we illustrate the adiabatic approximation on 
our branching-annihilation example. Here the peri- 
odic solution of the mean-field equation (fTT|) is q(t) ~ 
fl (1 + e cos uit): the mean-field fixed point follows the 
slowly changing rate A(t) adiabatically. 

The extinction rate of the branching-annihilation sys- 
tem with time-independent rates, was found, including 
the pre-exponent, in Refs. (l3l. \l4. [Til] : 



Q3/2 



(56) 



where S Q = 20(1 - In 2) and Q = A//i > 1. We intro- 
duce a slow sinusoidal modulation in the branching rate 
A, see Eq. (fT3)) , while keeping fx = //o constant. The 
instantaneous extinction rate becomes 



r ex (t) 



O 3 / 2 
2^ 



(1 +ecos^) 3/2 e" 5o(1+ecosii ' t) . (57) 



For this result to be valid it is necessary that the ar- 
gument in the exponent be large: <Sb(l + £cosujt) ^> 1, 
which yields 



1 - |e| > i/n. 



(58) 



The instantaneous extinction rate (|5T[) can be also ob- 
tained, with exponential accuracy, in the following way. 
One treats the time in the time-dependent eikonal Hamil- 
tonian (|14p as a parameter and looks for zero-energy tra- 
jectories as in the time-independent case. In this way one 
obtains the "adiabatic instanton" 



q = q a (p;t) 



2n(i+p) 

(2 + P) 



(1 + ecosuit) 



(59) 



and finds the adiabatic action So(t) = — J Q 1 q a {p\ t) dp — 
Sq(1 + s cos ut) which yields (the minus of) the logarithm 
of the instantaneous extinction rate, in agreement with 
Eq. (|57[) . This simple derivation yields a necessary con- 
dition for validity of the adiabatic approximation. In- 
deed, the adiabatic picture demands that u> be much 
less, at all times, than the adiabatically varying relax- 
ation rate A(i) = Ao(l + ecoswt). Therefore, we must 
require lj Aq(1 — |e|), that is 



a < 1 



(60) 



For small or moderate modulation amplitudes, \s\ < 1, 
one recovers the same low-frequency criterion a<l that 
led to Eq. (|3"4"|) . As |e| closely approaches 1, however, 
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FIG. 8: (Color online.) Comparison of the adiabatic instan- 
taneous extinction rate, given by Eq. (I57|l (the solid line) 
and the instantaneous extinction rate determined from a nu- 
merical solution of the master equation Eq. ©. The nu- 
merical extinction rate (the dashed line) was computed as 
-(d/dt)ln[l - Vo{t)]. The parameters are Q = 25, e = 0.2, 
and a = 0.04. 



the criterion (|60[) becomes much more stringent. The 
criteria ([55]) and ([50]) can be rewritten as 



1 - \e\ > max(fT 



(61) 



To verify Eq. (|57p. we solved numerically (a truncated 
version of) the original master equation |2|). The nu- 
merical instantaneous extinction rate was computed from 
-{d/dt) ln[l - V (t)]. Figured compares Eq. ([57]) with a 
numerical result for a small rescaled frequency a = 0.04 
and £ = 0.2, and good agreement is observed. The de- 
viation in the peaks is mainly caused by 1/ fi-corrections 
to the eikonal theory: we checked that it goes down as Q 
increases. 

The average extinction rate (|55[) can be written as 



47r 3/2 



2tt 



(1 + ecosT) 3 / 2 e _,So(1+ecosT) dr . (62) 



This expression is w-independent. For a comparison 
with the LT prediction (see below), it will be conve- 
nient to compute — ln(f ex ): the total action which in- 
cludes, in addition to the geometrical-optics contribu- 
tion, the physical-optics contribution coming from the 
pre-exponent. 

The integral in Eq. (|6"2"|) can be calculated analytically 
in two limits. For \e\D, 3> 1 one can employ the saddle 
point approximation and obtain 



n(i-M) 3/a 

47rVH(l-ln2)' 



-So(i-N) 



(63) 



The saddle point is located at t» = ±tt/uj (for e ^ 0, 
respectively); the effective width of the Gaussian is a ~ 
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FIG. 9: (Color online.) (a) Comparison between the action S 
given by Eq. (|64[) for eQ, S> 1 (the dashed line) and the natural 
logarithm of the MTE found by solving the master equation 
numerically (the solid line), as a function of e. (b) A blowup 
of the small-e region: shown are the eQ. <C 1 asymptote of 
the action S [Eq. (|66[) ] (the dash-dotted line), the eQ. S> 1 
asymptote of the action S [ Eq. (|64|l ] (the dashed line), and 
the natural logarithm of the MTE found by solving the master 
equation numerically (the solid line), as a function of s. The 
parameters in (a) and (b) are = 100 and a = 0.05. One can 
see that Eq. (|66|l breaks down when eQ, ceases to be small, 
whereas for eSl 1 the e-dependence becomes approximately 
linear, and the action is well described by Eq. (|64[) . 



(2|e|£lo; 2 ) 1 / 2 <C tt/lu. The resulting action is 



0.(1 - lei) 3 / 2 
5 -ln(f e:c )=So(l-kl)+ln - ) £U 

47r A /|e|(l - In 2) 

(64) 

The leading term Sq(1 — |e|) coincides with the zero- 
frequency limit predicted by the LT of Section IIIII see 
Eq. (|34|) with the a 2 -term neglected. Its physical mean- 
ing is transparent: in view of the adiabatically slow rate 
modulation, the effective "activation barrier" to extinc- 
tion iSo(l — |e|) is determined by the minimal value of 
A(t) = Ao(l + ecosut) which is equal to A (l — |e|). 
Equation (|64|) . however, also includes an important pre- 
exponent (recall that 17 3> 1), missed by the geometrical- 
optics order of the eikonal expansion. Note that, because 
of the additional pre-exponent, coming from the saddle- 
point integration, the ^-dependence of the resulting pre- 
exponent changes and becomes ~ f2 instead of ~ f2 3 / 2 . 
Fig. [9^ compares the prediction of Eq. (|64|) with the nat- 
ural logarithm of the MTE obtained by numerically solv- 
ing the master equation. 

In the opposite limit, \e\Sl <§C 1, the integral (|6"2")) can be 
calculated via a Taylor expansion of the integrand in efl. 
Terms proportional to cost vanish after the integration, 



12 



and the leading-order result is 



„(o) 



[l + e 2 Sl 2 (l -ln2) 2 ] 



This yields 



S = So 



e 2 n 2 (l 



In 2) 



(65) 



(66) 



Note that the strong inequality |e|0 1 in Eq. (|64|) 
coincides with the validity criterion of the eikonal method 
in the linear theory (LT): the left inequality in Eq. (|2"6"|) . 
Correspondingly, when the opposite strong inequality 
|e|f2 <C 1 holds, the correction to So that appears in 
Eq. (|66[) is much smaller than 1 and therefore strongly 
non-eikonal. Fig. shows a comparison between 
Eq. (f6"6"|) and the numerical solution of the master equa- 
tion, and excellent agreement is observed. 

One can notice that the leading-order term in the ac- 
tion (f6"l)) . 5o(l — |e|), goes down linearly in |e| for all 
\e\ < 1. The linearity in |e| is non-generic and stems from 
the fact that A(i) / /j,q is linear in e, and so is the exponent 
in the instantaneous extinction rate (|62p . In the general 
case, the |£|-dependence of the exponent that appears in 
the instantaneous extinction rate is nonlinear, and there- 
fore the linear ^-dependence of the action (which also 
appears in the LT-regime of the eikonal approximation) 
will break down at |e| = 0(1). This is indeed what one 
observes if, instead of modulating A, he modulates \x: 



X(t) = Xq , /i(i) = (1 + e cos cot) 



(67) 



In this case, using Eq. (|5"5")1 , we obtain the following action 
for \e\fl > 1: 



S = -ln(r ex ) = 



S 



1 



-In 



47rVN(l + N)(l-ln2) 

(68) 

For |e| <C 1 the leading-order term S /(l + |e|) ~ S (l — 
\e\) coincides with the leading term of Eq. or with 

Eq. where one must put a — 0. For not small |e|, 
however, the |e|-dependence is non-linear. 



V. SUMMARY AND DISCUSSION 

We have developed three complementary perturbation 
techniques for an analytical calculation of the exponen- 
tial reduction in the mean time to extinction (MTE) in 
single-species birth-death processes (or reaction kinetics) 
that occur in a time-periodic environment. These are the 
linear theory (valid at small modulation amplitudes), the 
Kapitsa method (valid at high modulation frequencies), 
and the adiabatic approximation (valid at low modula- 
tion frequencies). We presented our theory on the exam- 
ple of a simple branching-annihilation process. Figure [TU1 
shows different regimes of this process on the parameter 
plane (e > 0, a). 

The linear theory and the Kapitsa method, as we used 
them, are rooted in the time-dependent eikonal theory 
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t h e o r y '""„ 
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FIG. 10: (Color online.) A schematic diagram of different 
regimes of extinction on the parameter plane (e > 0, a) of the 
branching-annihilation process. Either the LT, or the Kapitsa 
method is applicable well to the left of and above the dash- 
dotted line (where the inequality A<S/5o <C 1 holds) and well 
to the right of the dashed line (where AS ^> 1). On the solid 
line e = /i(a)//2(a), see Eqs. (|47)l - (|49)l . The adiabatic ap- 
proximation holds well below the dotted line. The horizontal 
size of the non-eikonal regions at a — » is of order I/O. The 
rest of the diagram is self-explanatory. 



that employs a large parameter f2: the average number 
of individuals in the quasi-stationary state at times short 
compared with the MTE. The adiabatic approximation 
also employs the large parameter f2. It yields, however, 
a more accurate result than the other two techniques as 
it yields a nontrivial ^-dependent pre-exponent of the 
MTE. The higher accuracy, achieved in this way, is im- 
portant in view of the fact that, at low frequencies, the 
reduction in the MTE is the largest. 

An important result of this work is that, at high mod- 
ulation frequencies and not too small modulation ampli- 
tudes, the linear theory can greatly underestimate the 
true reduction in the MTE. In this regime an accurate 
result is given by the Kapitsa correction, calculated in 
this work. 

Our calculations of the modulation-induced reduction 
in the MTE are closely related to finding the optimal 
path to extinction, or instanton connection [T3, fll| . 
Although in generic Hamiltonians the existence of an in- 
stanton connection (a more standard mathematical term 
is a heteroclinic trajectory) can be g uaranteed only at 
very small modulation amplitudes [19( , our numerical re- 
sults strongly indicate that, in the birth-death processes 
with an absorbing state at zero, a perturbed instanton 
persists at any reasonable modulation amplitude and 
at any frequency. This remarkable fact must be inti- 
mately related to the non-generic structure of the birth- 
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death Hamiltonians. More specifically, a future proof of 
the instanton persistence at arbitrary (physically reason- 
able) time-dependent perturbations of the reaction rates 
will most likely exploit the invariance of the zero-energy 
mean-field line p — and the zero-energy extinction line 
q = in the perturbed Hamiltonian. 
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